Collective excitations in electron-hole bilayers 
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We report a combined analytic and Molecular Dynamics analysis of the collective mode spectrum 
of an electron-hole (bipolar) bilayer in the strong coupling quasi-classical limit. A robust, isotropic 
energy gap is identified in the out-of-phase spectra, generated by the combined effect of correlations 
and of the excitation of the bound dipoles; the in-phase spectra exhibit a correlation governed 
acoustic dispersion for the longitudinal and transverse modes. Strong nonlinear generation of higher 
harmonics of the fundamental dipole oscillation frequency and the transfer of harmonics between 
different modes is observed. The mode dispersions in the liquid state are compared with the phonon 
spectrum in the crystalline solid phase, reinforcing a coherent physical picture. 
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Electron/hole bilayer systems consisting of a pair of 
two dimensional (2D) layers of electrons and holes of 
equal density, (bipolar bilayers -BPBL-s) have attracted 
a great deal of attention over the past decade [l]. The 
main feature of the system is that in the appropriate pa- 
rameter domain, the electrons and the holes bind to each 
other in a dipole- like excitonic formation LlLyl- These ex- 
citons may also form a Bose condensate QMS 111- It is 
also expected that at strong enough coupling the system 
undergoes a transition into a Wigner crystal-like solid 
phase [2,13 ■ A recent work has also pointed out the excit- 
ing possibility of the existence of a supersolid phase over 
the dipole-solid domain |(|. A study of a strongly cou- 
pled classical BPBL, based on Monte Carlo simulation 
revealed the existence of dipole-liquid, dipole-solid, 
Coulomb-liquid and Coulomb-solid "phases", depending 
on the value of the coupling parameter and the layer sep- 
aration d. While the details of the boundaries and tran- 
sitions between the different phases would be different in 
the classical and quantum domains, there is little doubt 
that the topology of the classical phase diagram is quite 
generally correct Q. 

The problem we address in this Letter is the collective 
mode spectrum of the BPBL in the strong coupling (SC) 
regime, in a parameter range spanning all four phases 
but focusing on the more important and more intrigu- 
ing SC liquid state. The question of affinity of this spec- 
trum with that of the electronic bilayer (EBL), where the 
charges in the two layers are identical.arises; this latter 
system is now fairly well understood 0, H, El ■ From the 
point of view of the collective mode spectrum, the domi- 
nant feature of both systems in the SC regime is the local- 
ization or quasilocalization of the particles, either in the 
crystalline solid or in the SC liquid phase 10). This fact 
and the physical separation of the two oppositely charged 
layers that prevents their collapse, allows one to repre- 
sent the dynamics of the collective modes in the BPBL, 
even in the quantum domain, through a classical mod- 



eling The in-layer exchange effects, precisely 

because of the quasilocalization of the charges, contribute 
only to generating an exchange-correlation energy, which 
is well emulated by the equivalent classical correlation 
energy [l2T |. and tunneling is insignificant, because for 
any reasonable layer separation d value, d/as S> 1. This 
latter condition also ensures that the excitation energy 
of the exciton is well described in terms of its classical 
Kepler frequency. 

Our approach is based on a combined classical Molec- 
ular Dynamics (MD) simulation and theoretical analysis. 
To study collective excitations, simulation of a large as- 
sembly of particles is needed, which, at present, becomes 
feasible only by using a classical approach. What the 
classical approach obviously misses, however, is the effect 
of condensate on the dispersion. This is further discussed 
in our conclusions. 

While the hole and electron masses are, generally 
speaking, unequal, a simplified symmetric model with 
equal masses gives insight into the dominating aspects 
of the mode spectrum. Within the classical modeling, 
the symmetric BPBL system is then characterized by (i) 
the in-layer Coulomb coupling coefficient Y — (e 2 /akT) 
[a = {mr)^ 1 / 2 is the Wigner-Seitz (WS) radius, n is the 
areal density of particles] , and (ii) the separation of the 
two layers d. Addressing first the theoretical description 
of the collective mode spectrum, one should realize that, 
in contrast to the case of the EBL, a perturbative ap- 
proach is bound to fail - even at low Y values - due to the 
dominance of the excitonic bound state in the spectrum. 
Nevertheless, as a matter of orientation, we can start with 
the Random Phase Approximation (RPA) for the mode 
spectrum. The interaction between the particles can be 
characterized by the matrix (pn(k) = tp 2 2 (k) = 2ne 2 /k, 
<y9i 2 (k) = (f2i{k) = — (2ne 2 /k) exp(— kd). Its diago- 
nalization in the layer space leads to the identification 
of 2 longitudinal modes as in-phase and out-of-phase 
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(±) modes: = u>o\/ka [1 =p exp(— kd)] 1 , where 
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uq = \/2e 2 /ma 3 . The small-/c behavior of the + mode 
is acoustic, here with a slope s+ = uiavad, while the — 
mode has the typical 2D u> oc \fka~ behavior. Thus, the 
RPA misses, as it must, the appearance of the frequency 
of the intrinsic oscillation (excitation) of the dipole in 
the spectrum. In addition, the RPA it also unable to 
account for the SC collective effect, in particular, the 
generation of transverse shear modes. In an appropriate 
SC description, both of these effects should be correctly 
represented. The approach we adopt is the Quasilocal- 
ized Charge Approximation (QLCA) JhIEIj the method 
that we previously employed to identify the gapped ex- 
citation in the collective spectrum of the EBL p|- The 
QLCA approach has also been applied to a variety of 
Coulombic and other systems in the SC liquid state [14| , 
with results that have been successfully tested against 
the outcomes of computer simulations ja, E| and labora- 
tory experiments [T2I ITU Il6[ in both the classical and 
quantum domains. The longitudinal (L) and transverse 
(T) QLCA dielectric matrices in layer space Q 
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become functionals of the longitudinal (L) and trans- 
verse (T) projections of the dynamical matrix D^(k) = 

dlEq Q^fijio) [hij (|k-q|) -Sij^iViliQ^uiq)], 
hij(q) are the Fourier transforms of hij(r) — gij(r) — 1, 
where g%j{r) is the pair distribution function and A is a 
large but finite surface. The resulting mode structure 

c^(k) = [u%ka (1 T exp(-fcd)) + D ni k ) ± D U k )] V ' (2) 
^(k)=[^(k)±^ 2 (k)] 1/2 (3) 

now exhibits a behavior significantly different from its 
RPA counterpart. These dispersion curves are displayed 
in Fig^together with MD simulation results (see below). 

(i) At k = the — mode becomes a gapped mode (see 
also Fig.|2Ji), with 



co^(k = 0) = uj T {k = 0) = Q G (d) 



(4) 



where Qcid) is a functional of /ii2(r) as given in I n 
the BPBL hvz{r) is governed by a central peak h' 12 (r) 
around r = 0| (see Fig. Et). For d/a < 1 fl G (d) is 
well described by replacing h 12 (r) by h' 12 (r). When this 
central peak is approximated by a Gaussian (represent- 
ing a thermally excited dipole) fldd) becomes Q' K (d), 
the thermally broadened small amplitude oscillation Ke- 
pler frequency fluid) = \/2e 2 /d 3 m = LOo(a/d) 3 ^ 2 of the 
oscillating dipole. Dipole-dipole correlations shift flc(d) 
from this value only for d/a > 1. 

(ii) The + mode may be regarded as a density oscil- 
lation of the dipoles: with a 1/r 3 type dipole-dipole in- 
teraction, this suggests an w oc fc dispersion. This is 
borne out by the present calculation: the + mode shows 
an acoustic behavior, similarly to the RPA result, but 
with a phase velocity instead of being proportional to 
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FIG. 1: Dispersion relations for the L+, L-, T+, and T- 
modes (k = ka) at T — 40 and different layer separations. 
The legend shown in the d — d/a — 0.6 plot applies to all 
panels. Lines: QLCA, symbols: MD simulation. Inset shows 
the observed values of fccutoff of the T+ mode. 
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FIG. 2: (a) Dependence of the Kepler frequency Qk, the 
thermally shifted Kepler frequency Q' K and the gap frequency 
Qq on d at T = 80; also shown for comparison is the EBL 
gap. (b) s+ sound velocity vs. d for F = 80. (c) MD pair 
correlation functions for F — 80 and d = 0.6 (left and right 
vertical scales belong to the left and right horizontal parts). 



{d/a) 1/2 , is of the order of d/a: s = oj+(k — > 0)/k = 

uj d [(99/96) / dfg 12 {r)/f 2 ] 1/2 (see Fig.Eb)- 

(iii) In addition to the 2 longitudinal modes 2 trans- 
verse shear modes appear, with a behavior qualitatively 
similar to their longitudinal counterparts, except that, 
as expected for shear modes in a liquid, the T + mode do 
not extend below a finite fc cu toff value (see inset in Fig. 

m 

(iv) For k — ► oc all four modes approach the Einstein 
frequency (the frequency of oscillation of a single particle 
in the frozen environment of the others) [?J of the system: 
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where £1% 2D 



(wg/2) / dfgnir)/f 2 = 0.38w§ is the Ein- 
stein frequency of an isolated 2D layer 01 ( see later). 

Illustrative longitudinal and transverse current-current 
correlation spectra obtained from our MD simulations 
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FIG. 3: Longitudinal current (a,b) and transverse current 
(c,d) fluctuation spectra obtained at Y — 40 and d = 0.6. 
The arrows indicate increasing wave numbers k as listed in the 
panels. Note the appearance of higher harmonics of the gap 
frequency: the second and fourth harmonic in the + spectra, 
and the third harmonic in the — spectra. 



fundamental) in the + spectra) [Fig. ^ . These harmonics 
are the most pronounced around T = 40 in the strongly 
coupled dipole liquid phase, showing a diminishing trend 
both for lower and higher coupling values: at high T 
values, the low amplitude of the thermally excited oscil- 
lations is not conducive to the generation of harmonics, 
while at lower T values, the thermal motions damp the 
higher harmonics. The analysis of the details of the non- 
linear processes underlying the harmonic generation will 
be presented elsewhere. 

We calculate the phonon spectrum of the crystalline 
phase, through the standard harmonic approximation, 
by summing over a lattice of 2 x 10 7 sites. In the solid 
phase the BPBL crystallizes in a hexagonal structure, 
with the particles in the two layers facing each other (for 
any layer separation, in contrast to the EBL which ex- 
hibits a variety of lattice structures) . One can recognize 
the 4 modes identified in the liquid state, with an ex- 
pected anisotropic dispersion and with the understanding 
that the "longitudinal" and "transverse" labeling of the 
modes describes their polarization for propagation along 
the principal crystal axes only. 



(based on the Particle-Particle Particle-Mesh (PPPM) 
method with periodic boundary conditions) are dis- 
played in Fig. for T = 40 and d/a = 0.6. Collec- 
tive modes show up as peaks in the correlation spectra. 
The MD-generated dispersion relations inferred from the 
spectral peaks for the 4 modes for T = 40 and for a se- 
ries of d/a values are displayed in Fig. ^ together with 
theoretically calculated values. Agreement between the 
QLCA theory and simulation data is very good, espe- 
cially for lower k values. Deviation from the character- 
istic <i 3 / 2 behavior of £Ik{<1) seems to occur on passing 
the (dipole liquid) /(Coulomb liquid) phase boundary Q- 
Comparison with the EBL spectrum reveals a more 
pronounced and more robust gap value (Fig. Eh) i which is 
not unexpected if one is mindful of the different physical 
mechanisms responsible for creating the gap. 
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FIG. 4: Longitudinal current fluctuation spectra showing the 
higher harmonics of the gap frequency for k — 3.54, d — 0.6, 
r = 10, 40 and 120. The transverse spectra exhibit similar 
behavior. 



A remarkable feature of the simulation spectra is the 
appearance of higher harmonics of the gap frequency: we 
observe the emergence (a) of the third harmonic in the 
— spectra and (b) of the second harmonic (but not the 
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FIG. 5: (a) Lattice dispersion of the quasi-L+ and quasi-L_ 
(full lines) and quasi- T + and quasi- T_ (dashed lines) modes 
for a series of lattice angles between 0° < a < 30° (a = 
0° corresponds to the nearest neighbor direction, the arrows 
indicate increasing a), (b) Comparison at d — 0.6 of QLCA 
(at r = 80) and angular averaged lattice dispersion. 

The mode spectrum is shown in Fig.03 For the sake of 
comparison with the QLCA dispersion, also shown are 
the angle-averaged longitudinal and transverse projec- 
tions: the similarity between the two is the consequence 
of the underlying microscopic order in the SC liquid as 
predicted by the QLCA. The — modes in the lattice spec- 
trum exhibit a k = gap, very similar to the one in the 
liquid spectrum; in contrast to the EBL, the gap is not 
split by the lattice anisotropy. This is a consequence of 
the quasi-isotropy of the hexagonal lattice environment 
to 0(k 2 ). Physically, the — modes can be regarded, in 
a good approximation, especially for small d/a values, 
as the superposition of a + mode on the Qk dipole os- 
cillation frequency. The angle-averaged modes emulate 
the liquid-phase mode dispersion for low and moderate 
/c-values, but deviate substantially from it for k — ► oo 
(see Fig. Eb). To understand this difference, one has 
to focus on h' 12 (r), which shows a quasi-Gaussian be- 
havior in the liquid, but becomes a delta-function in a 
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perfect T = lattice. The contribution of oscillating 
dipoles within a Gaussian distribution results in a su- 
perposition of random phases at k — > oo, averaging out 
to 0, ensuring the limit shown in Eq. [S] In contrast, the 
delta-function contribution generates additional coherent 
k- independent terms that survive for k — > oo. The simu- 
lation results clearly indicate that even for high T values 
the large but finite k dispersion curves follow the trend 
set by Eq. describing the liquid, rather than the one 
corresponding to the lattice. The great sensitivity of the 
behavior of the mode dispersions at higher k values to 
the width of h' 12 (r) could motivate an observational tech- 
nique for acquiring information of the width of central 
peak distribution through tracking the high-/c behavior 
of the collective excitations. 

The phonon spectrum in the BPBL was the subject of 
a recent investigation by Ref . [2lJ . This work addressed 
the behavior of the dipolar "phase" , within an approxi- 
mation geared to the assumed separability of the overall 
collective and intrinsic dipolar excitations: our results 
show that this technique, reasonable as it may appear, 
leads to an unphysical behavior in the calculated mode 
dispersion. In particular, pH finds (i) that the gap is 
split and there are distinct transverse and longitudinal 
gap values; (ii) that both of these exhibit an anisotropic 
behavior; (hi) that for k — > the — mode dispersion 
curves exhibit a negative slope; and (iv) that the k — > oo 
behavior does not conform to that demanded either by 
Eq. 0or by the lattice condition. These features violate 
basic physical principles pointed out above. 

In summary, we have determined the dispersion char- 
acteristics of a strongly coupled symmetric bipolar (elec- 
tron/hole) bilayer in the strongly coupled domain, where 
the dominant exchange-correlation energy provides a suf- 
ficient quasilocalization of the particles to engender a 
quasi-classical behavior. Our analysis in the strongly 
coupled liquid phase is based on the MD simulation of the 



density and current fluctuation spectra and on the appli- 
cation of the theoretical QLCA technique. The bench- 
mark phonon spectrum in the crystalline solid phase is 
determined through lattice summation technique. The 
results of all these various approaches reinforce each other 
and provide a coherent physical picture. The 4 modes, 
characteristic of bilayer systems (in-phase/out-of-phase, 
longitudinal/transverse) emerge, with the expected k = 
energy gap in the spectrum of the out-of-phase mode 
0, @j- The latter, is now mostly governed by the in- 
trinsic oscillation frequency of the dipoles and affected 
by collective interaction only in the Coulomb phase, for 
d > a. A remarkable effect, unique to the strongly cou- 
pled liquid phase, the generation of harmonics of the gap 
frequency and the transfer of the even harmonics from 
the out-of-phase mode to the in-phase mode, has been 
observed. 

The presence of the condensate would manifest itself 
primarily by reducing the quasi-localization of the parti- 
cles (as clearly demonstrated by 22] ). A rough estimate 
of the degree of this reduction is provided by replacing 
h(k) (and thus D{k)) by h{k) = (1 - f) 2 h(k\ where / 
is the condensate fraction, as suggested by [23. Since / 
in the SC regime is expected to be quite small (/ = 0.02 
@), this should not be a dramatic effect. In the solid 
phase, if the supersolid develops, the system would prob- 
ably behave as a binary mixture of the normal crystal 
lattice and the superfluid liquid constituted by the va- 
cancies |2 1[ , with a coupling between the lattice phonons 
and the collective excitations of the latter. These issues 
have to be addressed elsewhere. 
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